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Abstract 

Competing styles in Statistical Mechanics have been introduced to investigate physico-chemical 
systems displaying complex structures, when one faces difficulties to handle the standard formalism 
in the well established Boltzmann-Gibbs statistics. After a brief description of the question, we 
consider the particular case of Renyi statistics whose use is illustrated in a study of the question of 
the "anomalous" (non-Fickian) diffusion that it is involved in experiments of cyclic voltammetry 
in electro-physical chemistry. In them one is dealing with the fractal-like structure of the thin 
film morphology present in eletrodes in microbatteries which leads to fractional-power laws for 
describing voltammetry measurements and in the determination of the interphase width derived 
using atomic force microscopy. The fractional-powers associated to these quantities are related 
to each other and to the statistical fractal dimension, and can be expressed in terms of a power 
index, on which depends Renyi's statistical mechanics. It is clarified the important fact that 
this index, which is limited to a given interval, provides a measure of the microroughness of the 
electrode surface, and is related to the dynamics involved, the non-equilibrium thermodynamic 
state of the system, and to the experimental protocol. 

KEYWORDS: Complex structured systems, anomalous diffusion, fractal surfaces, unconventional 
statistics, Renyi's statistics. 
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I. INTRODUCTION 

A couple of decades ago Montroll and Shlesinger [l| stated that in the world of inves- 
tigation of complex phenomena that requires statistical modeling and interpretation, several 
competing styles have been emerging, each with its own champions. Lately a good amount of 
efforts have been devoted to the topic. What is at play consists in that in the study of certain 
physico-chemical systems we may face difficulties when handling situations involving fractal- 
like structures, correlations (spatial and temporal) with some type of scaling, turbulent and 
chaotic motion, small size (nanometer scale) systems with eventually a low number of degrees 
of freedom and complicate boundary conditions, and so on. The interest has been nowadays 
enhanced as a result that such situations are present in electronics, opto-electronics, etc, 
devices of the present day point-first technologies, on which is so dependent our society, 
also in disordered systems, polymeric solutions, ion conducting glasses, microbatteries, and 
others. Most cases involve, what in the literature on the subject are called "anomalous" 
physical properties ( let them be transport, optical, hydrodynamical, etc. ones). As a rule 
what is noticed is that the observable experimental data can be theoretically adjusted by 
means of fractional-power laws in place of the standard ones. We can mention one type, 
related to the one presented here, consisting into "anomalous" (non-Fickian) diffusion in 
micelles ^ 0, |^ , where the analysis of data used to be done introducing the so-called 
Levy statistics. In this paper we present a closed statistical theory of a large scope in place 
of such approach. 

In these situations the difficulties consist, as a rule, in that the researcher is impaired, for 
one reason or other, to satisfy Fisher's criteria of efficiency and sufficiency p in the standard 
formalism to implement the conventional, well established, physically and logically sound 
Boltzmann-Gibbs statistics, meaning an impairement to correctly include the presence of 
large fluctuations (and eventually higher-order variances) and to account for the relevant and 
proper characteristics of the system, respectively. To contour these difficulties, and to be able 
to make predictions - providing and understanding, even partial, of the physico-chemical 
properties of the system of interest, for example, in analyzing technological characteristic 
of a device - there has been developed efforts along the line of introducing heterotypical 
statistics for providing an auxiliary method for fulfilling such purposes. 

Among existing approaches it can be mentioned, what we call Generalized Statistical Me- 
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chanics as used by P.T. Landsberg showing that functional properties of the (informational) 
entropy (which is a generating functional of probability distributions) provides different 
types of thermodynamics, and rise the question of how to select a "proper" one, that is, 
are some better than others? 3]; Superstatistics as used by C. Beck and E.G.D. Cohen for 
nonequilibrium systems with complex dynamics in stationary states with large fluctuations 
on long-time scalesla, y. 

n 

So-called nonextensive statistics based on Havrda-Charvat statistics [10|| applied to a 
number of cases are described in the Conference Proceedings of Ref. [111. Renyi statistics 

nn n 

12, 1^1 has been introduced in the scientific literature [1^, with T. Arimitzu and P. Jizba 
ISL llq recently introducing an extensive analysis in a paper intitled " The world according 
to Renyi; Kappa (sometimes called Deformational) statistics - a particular case of Sharma- 



Mittal statistics 



17| - was used in problems of plasma physics in celestial mechanics |l8| 



and by G. Kanadiakis in the case of relativistic systems | 

The introduction of these alternative statistics resides on the fact, already stated, that the 
well established mechanical-statistical treatment via Boltzmann-Gibbs ensemble formalism, 
let it be around equilibrium (as in the usual response function theory) or for systems arbi- 
trarily away from equilibrium [2^ , has its application impaired because of the complications 
set forth by the presence of spatial correlations, resulting for example in the illustration 
here presented from spatially varying microroughened boundaries showing a fractal struc- 
ture. Such fractal structure introduces a degree of fuzziness in data and information, which 
consists in the lack of information on the space correlations imposed by the fractal-like mi- 
croroughness of the electrode surface. Hence, as noticed, it is not satisfied Fisher's criteria of 
efficiency and sufficiency in the characterization of the system, that is, the chosen description 
does not summarize the whole relevant information for the situation in hands. Hence the 
use of formulations based on unconventional statistics 22] to perform calculations of average 
values and response functions in which the accompanying heterotypical probability distri- 
butions preferentially weight the probabilities of the states that contribute to the dynamics 
involved We restate that the use of unconventional (index dependent) heterotypi- 

cal statistical mechanics for systems displaying correlation effects of one kind or other, is a 
circumventing way for overcoming the difficulty to assess such characteristics that should be 
incorporated in the conventional formulation within the well established Boltzmann-Gibbs 
theory. 
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We consider here Renyi statistics [12, ll3| applied in statistical mechanics [15|, |16|, |23|, |2J| , 
illustrating its use and utility in the case of a particular experiment, namely, observation 
and measurements using cyclic voltammentry in the study of fractal-structure electrodes 
in microbatteries. First, we begin again noticing that the advanced technologies that keep 
evolving nowadays present challenges for the associated physics which describes the phenom- 
ena that are in them involved thus, requiring a careful use of the accompanying theoretical 
physics to face the difficulties associated with the handling of these new situations. They 
involve aspects of the physics of soft-matter and disordered systems, complex structured sys- 
tems, the already mentioned small systems, and others, as well as the question of functioning 
in far-from-equilibrium conditions and in the presence of ultrafast (pico- and femto-second 
time scales) relaxation and pumping effects. We are dealing in all these cases with condensed 
matter and then, we recall, it is inevitable for their study at the physico-chemical level to 
resort to statistical mechanics - in a sense it can be said that the physics of condensed matter 
is statistical mechanics by antonomasia. We consider here a particular situation which pro- 
vides an excellent illustration of what can be called Unconventional Statistical Mechanics at 
work. We begin noticing that as a consequence of the nowadays large interest associated to 
the technological applications in displays, electrochromic windows and microbatteries (see, 
for example, Ref. [25]), the study of growth, annealing and surface morphology of thin film 
depositions has acquired particular relevance. These kind of systems involve the presence 
of microroughened surface boundaries in a geometrically constrained region, when fractal 
characteristics can be expected to greatly infiuence the physical properties ^26]. It can be 
recalled that fractals display self-similarity, but the latter is said to come in two flavors: 
exact and statistical, the latter being the one with which we are dealing here. Moreover, 
fractals are characterized in terms of their, in general, fractional "dimension" - as we shall 
see in the present case - which is a "measure" of the complexity of the structure. 

Hence, we do have that the dynamics or hydrodynamics involved in the functioning of 
such devices can be expected to be governed by some type or other of scaling laws j^, 2^ ■ 
Such characteristics have been experimentally evidenced and the scaling laws determined 
by researchers who resorted to the use of several experimental techniques, as atomic force 
microscopy 30j which allows to obtain the detailed topography of the surface - and then 



measuring the statistical fractal dimension - and cyclic voltammetry 



which is an elec- 



trochemical technique used for the study of several phenomena, also allowing for the charac- 
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terization of the fractal properties (scaling laws) of the system. Particularly, through cyclic 
voltammetry it can be evidenced the property of the so-called "anomalous" (more properly 
called non-Fickian) diffusion in these systems. It is our aim to present here a theory appro- 
priate to deal with these kind of situations. For dealing with transport processes in elec- 
trochemical cells we introduce two different aspects of recent thermo-hydrodynamic theory: 



one is Informational Statistical Thermodynamics (1ST for short 1321 l33l |. which provides sta- 



tistical foundations to Extended Irreversible Thermodynamics |3J|), and a non-equilibrium 



ensemble formalism [2l|. However, for the kind of experiments to be analyzed one needs, 
as noticed, to resort to a modified form of the standard formalism in the original and well 
established Boltzmann-Gibbs foundations of statistical mechanics. Originally, the relation 
between irreversible thermodynamics and transport laws, as Fick and Fourier laws, were well 
established - and kinetic and mechanical statistical foundations given - in the framework 
of classical (or Onsagerian) irreversible thermodynamics. But, it faces limitations since the 
theory is restricted to weak amplitudes of movement, with a linear relationship between 
fluxes (currents) and thermodynamic forces, smooth variation in space and time, locality in 
space and instantaneously in time (i.e. neglecting space and time correlations) and weak 
fluctuations, and then does not cover the case of fractal-structured systems. To remove 
these restrictions appeared a variety of approaches, a most convenient one is the already 
mentioned Extended Irreversible Thermodynamics and its statistical counterpart Informa- 
tional Statistical Thermodynamics, and for a practical handling of fractal and nanometric 
systems one can resort to unconventional statistics. 1ST is founded on a non-equilibrium 
ensemb 
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e formalism which provides a generalized nonequilibrium grand-canonical ensemble 
3^. On the other hand, unconventional statistics applies to systems displaying long- 
s'. It has been 



range interactions, persistent memory, evolution in a fractal space, etc [23, 
applied to some kind of studies of "anomalous" diffusion, as for example in Refs. [37-39], and 
we reconsider the question here, specified to the situation in fractal-structured electrodes of 
microbatteries leading to "anomalous" results in experiments involving cyclic voltammetry. 

II. THEORETICAL ANALYSIS 

Within the variational approach to Statistical Mechanics jj21j the use of unconventional 
(index-dependent) informational entropies (better called measure of uncertainty of informa- 
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tion, and not to be confused with the physical entropy which is a property of the thermody- 
namic state of physical systems, while the informational one is a property of any probability 
distribution j2^;Q;l3) 13 leads to the construction of alternative Statistical Mechanics 
as described in Refs. [23,24]. Basically t he p rocess consists of two steps: (l)The choice of 



probability 




an heterotypical probability distribution |22l l23l |2J|, say Qh, and (2) the use of an escort 



41[ of order 7 in terms of this heterotypical distribution, say V.y{gh}, namely, in 



the calculation of average values it is used the definition 

/^A)=Tr{Av,{g,}] , (1) 
where the escort probability of order 7 is given by 

where 7 is a real positive number (see for example Ref.[41], where the name escort probability 
is used, and the concept generalized, and also in Chapter IX, pp. 569 et seq. of Ref.[13]). 
An in depth presentation and discussion is given elsewhere 13 ' ^^^^ suffices to say 
that what the heterotypical distribution does is providing a weighting of the probabilities of 
the states involved, and the escort probability accounts for the influence of correlations and 
higher-order variances. 

Acccd,,, to BecU and Q , the .se of tKe esco. p.obab.Hty ..t.oduces an .nc.ease 

in information (by incorporating self-consistently correlations and variances) in the sense 
that, for the case of index a-dependent Renyi distribution (see below) 

(1 - «)' ^ = Tr {V^ {g^} [In {V^ {g^}) - In (g^)]} , (3) 

where la is the quantity of information (the negative of Renyi's informational entropy in 
terms of ga), and the right-hand side is interpreted as the information gain when using the 
escort probability built in terms of the original one g^- 

Let us now proceed to deal with the problem on hands, when we make use of an unconven- 
tional statistics in Renyi approach in order to provide a statistical-thermodynamic treatment 
of the "anomalous" diffusion law proposed to account for the experimental results obtained 
using cyclic voltammetry in the system of fractal-structured electrodes in microbatteries. It 
involves the process of motion of charges in the electrolyte bordered by a fractal-structured 
thin film on the electrode, motion which is usually considered that can be described as a 



process of diffusion governed by Fick law 31]. But the use of standard Fick diffusion equa- 
tion does not properly describe the experimental results, which show some kind of power 
law (see below) differing from the expected one. This is a consequence that in a constrained 
(nanometric scaled) fractal-like geometry, Fickian-diffusion does not apply because are not 
satisfied the quite restrictive conditions mentioned in the Introduction. A proper description 
using the conventional approach in statistical mechanics for providing foundations for hydro- 



dynamics requires the use of a generalized higher-order thermo-hydrodynamics [42 
that is, the description including as basic macrovariables the density of energy, the den- 
sity of particles, h{T,t) and n{r,t), their fluxes of all orders, Ih{r,t), In{r,t), ln\r,t) with 
r = 2,3,... indicating the order of the flux and its tensorial rank; all the quantities are 
defined in Refs.[21] and [42]. To resort to this higher order hydrodynamics (which provides 
at the microscopic, i.e., statistical mechanical level, generalizations of super-Burnett and 
super- Gruyer-Krumhansl equations (see Refs.[45] and [38]), even in the linear aproximation 
and neglecting fluctuations is a formidable task in the present case with spatially varying 
boundary conditions (in a nanometric scale and with a complex structure). Moreover, even 
;houKh a contraction in the choice of the basic set of macrovariables can be introduced 
2l[ [45^ , it is difficult to establish the order of the truncation. Therefore, what can be at- 
tempted in order to obtain a description of some properties of the system, is to introduce 
an unconventional statistics, leading to an unconventional hydrodynamics, in terms of a 
low-order truncation, in which one keeps only the variables 

{/i(r,t),n(r,t),I„(r,t)}, (4) 

that is, the energy density, the density of particles and the first-order flux (current). In that 
way we are led to obtain an "anomalous" diffusion equation depending on a fractional-power 
which is determined by the dynamics, the fractality, the geometry and dimensions, and the 



le unconventional generalized higher-order thermo- 
i^Ql)- Let us proceed to derive their equations of 



thermodynamic state of the system (t 
hydrodynamics is presented elsewhere 

evolution using Renyi approach to unconventional statistics |23l 24 1. For that purpose, first 
we separate the Hamiltonian in the form 

H = Ho + H', (5) 

where Hq is the kinetic energy and H' contains all the interactions that are present in the 
system. Moreover, in a classical approach, the basic dynamic variables, whose average values 



are those of Eq. (jH) , are 



Mr|r)= / d'p^n,ir,p\T), (6) 



2m 



n(r|r) = J d'^pni{r,p\T) , (7) 

Ur\T) = J rf3pZ^,(r,p|r), (8) 

where m is the mass of the particles, p the momentum and fii is the reduced single-particle 
density function 

N 

rii(r,p|r) = 5^5(r-r,)5(p-p,), (9) 

j=i 

with N being the number of particles, and F a point in phase space. The corresponding 
macrovariables are the average of those of Eqs.® to (jHI) using the indexed Renyi statistical 
approach, that is 

h{r,t) = J dTh{r\T)V^{T,t), (10) 
n{r,t) = J dTn{r\T)V^{T,t), (11) 

4(r,t) = j dr4(r|r)p,(r,t), (12) 

where the integration is over the phase space. Moreover, according to the theory ^ , 
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M 

2J| it is used the corresponding escort probability distribution (cf. Eq.(All) in Appendix 



A) 



V^iT,t) = g2iT\t) / I dVg'^^iTlt), (13) 



and it can be noticed that the order of the escort probability is the same of the index 
in Renyi' s distribution jl^, so we keep only one parameter in the theory (the so-called 
infoentropic index a), with Renyi's heterotypical distribution of probability ^)2(r|t) given in 
Eq.(A13) after using Eq. (A14), both in Appendix A. 
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The equations of evolution for the basic variables associated to the material motion (the 

nn 

one of interest for us here) derived in the context of a nonlinear kinetic theory |21l . |42| are 



^n(r,t) = -V-I„(r,t) 



|^I„(r,t) = -V-/f (r,t) + J„,(r,t) 



(14) 



(15) 



with Eq. p4|) being the conservation equation for the density, the terms with the presence of 
the divergence operator, V ■, arise out of the contribution resulting from performing the, in 
this classical case, Poisson bracket with the kinetic energy operator (i-e- {^^j-^o}, etc), 
and 

ft 



t) = \\mj dV dte<'-'^ [[Uv\^),H'{T\{t -t)] H'{T)] P„(r,t), (16) 



is a scattering integral accountin g fo r the effects of the collisions generated by H\ and /'^^ 
is the second-order flux given by 2l|, |36 

/r(r,t) = 



rf3plPPlni(r,p|r)D,(r,t), 



(17) 



where [...] stands for tensorial product of vectors, rendering a rank- two tensor. To solve the 
system of Eqs. (fT^ and (fT3j) we need in Eq. (fT3j) to express the right-hand side in terms of 
the basic variables. The scattering integral takes in general the form of a relaxation-time 
approach , namely (see Appendix B)) 



(r,t) 



InM) 



(18) 



where r is the momentum relaxation time (see for example Ref.[49]). Transforming Fourier 
in time Eq.lfT^ we have, after using Eq.lfTHj). that 



(1 + iuTio) In(r,t^) = -t/qV ■ /Ji^^(r,t^) , 

which in the limit of small frequency, meaning wt/q, -C 1, reduces to 

I„(r,t) = -r,,V-/^l(r,t), 

and then, after using Eq.ipUj) in Eq.lfT^. we obtain that 

d 



dt 



n 



(19) 



(20) 



(21) 



In order to close Eq. ()2H) we need to express the second-order flux in terms of the basic 
variables, n and I, which after some calculus, described in Appendix A, results in that 

V-V-/f (r,t) = UVV(r,t), (22) 

with C,na given in Eq. ()A20|) and where 7q, is the a-dependent fractional power 

5 — 3q! 

7a = ^ , 23 

6 — a 

and consequently it follows the so-called "anomalous" diffusion equation 

|^n(r,t) + D„„VV(r,t) =0, (24) 

where Dna = ^na^ia (it can be noticed that the dimension of D is cm^'''"~^/sec). Moreover, 
the possible values of a having physical meaning (i.e. not leading to singularities) belong to 
the interval [see Appendix A] 

l<a<^. (25) 

At this point we recall, and stress, that the "anomalous" diffusion equation is an equation 
appropriate for approximately describing the hydrodynamic motion, being an artifact (thus 
the use of the word anomalous between quotation marks) of having applied an unconven- 
tional statistics in the truncated description as given in Eq. (4). As already noticed, in the 
conventional (well established and parameter free) formalism for satisfying the principle of 
sufficiency it must be used a higher-order hydrodynamics. 

We proceed to analyze this question of "anomalous" results in measurements using cyclic 
voltammetry in fractal-structured electrodes in terms of the previous results. Let us, for 
the sake of completeness, briefly summarize the case: The difference of chemical potential 
between an anode and a cathode with a thin film (nanometric fractal surface) of, say, nickel 
oxides, produces a movement of charges in the electrolyte from the former to the latter. 
In a cyclic voltammetry experiment, these charges circulate as a result of the application 
of a potential V{t) , with particular characteristics: It keeps increasing linearly in time as 
Vq + vt, where w is a scanning velocity, during an interval, say At, and next decreases with 
the same sweep rate, i.e. as Vq + v A t — vt, until recovering the value Vq. A current i{t) 
is produced in the closed circuit, which following the potential V{t) keeps increasing up 
to a peak value ip , and next decreases. This current is the result of the movement of the 
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charges that keep arriving to the thin film fractal-hke cathode. It is found that there follows 
a power law relation between the peak value, ip, of the current and the rate of change v 
of the electric field, namely ip ~ w'^. The value of the current to circulate in the cathode, 
i{t), is proportional to the charge at the interface, that is, the value at it of n{x, t), to be 
given by the solution of Eq. (1261) which in one dimension, say, the x-direction normal to the 
electrode surface, is given bv Isoi Isij l 

n{x, t) = 6„ t-^'^ [a^ + t'^^'^] ^ , (26) 

where a (which ensures the normalization) and ba are constants of no specific interest in 
what follows and then we omit to write them down, and 

3 — a 
2 - a_ 

Hence the charge at each given point x at the interface is given by Eq.(26), and the 
current to be produced in the electrode - motion of this charge arriving at it under the 
action of the field to be applied - is proportional to it, and then we have the law 



/^a = [7a + 1] ^ = ^ 



(27) 



m = a t"''" [a' + t-'^^] - ;L./ace , (28) 

where is a constant. 

Taking into account that fia is positive (cf. Eqs. (27) and (25), and that < 7q, < 1), 
for not too-short times after application of the field we can expect that ^ x'^t^'^'^" and 
then there follows a power law in time for the current, namely 

i{t) ^ f-''" = = t-ii^ . (29) 

But taking into account that the applied potential is V = Vq + vt (where v is the scanning 
velocity), and then t = {V — Vo)/v, Eq. (29) leads to the power law 

ip^ v\ (30) 

where e is fia, that is, 

1 ^'^ - a' 

(31) 



e = (7a + 1)-^ = \ 



2 - a_ 

In the conditions of the work of Pajkossy and Nyicos '26], the fractal power e is related 
to the fractal dimension df of the electrode rough surface by the relation e ^ (d/ — 1)/2, and 
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then in experiments when it holds, from the log-log plot of ip vs v there follows the value of 
e and, consequently, the fractal dimension can be estimated. Using Eq. ()31|) we arrive to the 
relation between fractal dimension and the parameter a, namely a = (8e — 3)/ (4e — 1), and 
then 

(-) 

We can see that for df = 2 {a. perfectly flat surface) we have a = 1 as it should, while 
for the other possible extreme limit of df = 3, it follows a = 5/3. As already noticed, 
the theory precisely restricts the values of a to this interval; outside it the calculations 
present singularities. Then as a consequence the exponent e is also limited, and we can 
only expect for it values in the interval 0.5 < e < 1. Hence, we stress that these are the 
permited limiting values of the indexes a and e: the value which adjusts the experimental 
data is to be contained between these limits, being dependent on the systems dynamics, 
geometry and size, macroscopic thermodynamic state, and the experimental protocol (in 
the present case depends on the scanning velocity v, as shall be shown in the next section). 
We can see that the three quantities, a,e and df are related and give a kind of measure 
of the microroughness of the electrode surface. A perfect one (a totally smooth surface) 
corresponds to a = l,e = 0.5 and df = 2, and as the surface becomes more and more 
imperfect, these quantities increase, above those values, but with upper boundaries, being 
1.666 for a, 1 for e and 3 for df. 

On the other hand, atomic force microscopy allows to obtain the fractal dimension and 
a kind of degree of microroughness of the surface. According to the method {2^, it is 
calculated the so-called interface width, namely 



W{L) 



L 



(33) 



where L is the number of partitions in which the surface is divided and hi — his the deviation, 
within each partition, of the height with respect to the average value. There exists a power 
law between W and L, namely W ^ U , where r is called the roughness exponent, which 
is related to the fractal dimension as df ^ 3 — r. Therefore, comparing with the theory 
developed for the case of cyclic voltammetry, we obtain relationships relating the different 
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quantities, namely 



2(1 



5 — 3a 
2 - a 



27» 
7a + 1 



(34) 



Recalling that, according to the theory, 1 < a < 5/3, we have that r can only attain 
the values < r < 1, comprising the conditions between highest roughness to a perfectly 
smooth surface, and thus can be used as a kind of "order parameter", or, better to say, 
"degree of fractality parameter". 



III. THEORY AND EXPERIMENT 

I — I 

We apply the preceding theory to the study of a class of experiments [52j. In that work, 
the fractal dimension of nickel oxide thin films was analyzed by means of cyclic voltammetry 
and compared with the results obtained by atomic force microscopy (AFM). 

Figure 1 shows the dependence of the peak current with the sweep-rate velocity v. We 
can notice that the peak value of the current increases continuosly with and also that it 
can be identified two regions where the dependence is very approximately linear (one in the 
range 1 < f < 10 mVs"^ with a more pronounced slope that the other in the range 20 < u 
up to ~ 100 mVs~^). Once we do have a log-log plot, consequently it can be identified two 
values of the exponent e, in accordance with Eq.ijSOJ. Moreover, from equation Eq. (jHT|) we 
obtain that 

8e-3 

« = 4731 ' (35) 

and two values of the infoentropic index a can also be determined. The calculated values 
are presented in Table I. The associated statistical fractal dimension values (rf/, Eq.()32j)) 
are also presented in this table. The larger df values are obtained in the slow-sweep ranges, 
in comparison with those obtained at faster scans. The fact that can be identified from 
this is that the fractal dimension depends on the observational scale, which is limited by 
the experimental protocol, here the scanning velocity range: At low values of w, the ions 
are constrained in a spatial region close to the electrode (known as the diffusion layer), 
and are able to sense the fine details of the microroughness at the interface. On the other 
hand, increasing the scan rate also increases the length of the diffusion layer, and the global 
morphology of the electrode surface becomes predominant. This behavior is particularly 



13 



expected in the case of sensing by some experimental method the surface morphology of the 
thin film: this surface presents grains, and two distinct statistical processes are responsible 
either for the formation of the local surface morphology of individual grains or the global thin 
film surface morphology, formed by the assemblage of all grains. The experimental values 
determined in this work for the infoentropic index a are within the expected theoretical 
range (1 < a < 5/3), and the values obtained by cyclic voltammetry (Table 1) indicate 
that a stronger asperity can be associated to the local grain morphology, in constrast to a 
smoother thin film global morphology. 

Figure 2 presents a typical AFM micrography. The presence of structured grains, with 
mean size ?«180 nm, can be clearly seen. Figure 3 presents a log-log plot of the interface 
width WiV) as a function of the partitions L. Two scaling regions (at short and medium 
L ranges) are evidenced, separated by a transition L value of ^140 nm , close to the mean 
grain size. The roughness exponent r values calculated in accordance with Eq.()34|). and the 
associated AFM images statistical fractal dimension values are listed in Table I. The calcu- 
lated statistical fractal dimensions from both techniques are in good agreement. Returning 
to Fig.l, it is worth noticing that the full curve of dependence of ip with u, once we take 
into account Eqs. ()30|) and ()31|) . can be approximated by the law 

a[v) ~ , (36 

where v\ = 235 ± 12 and V2 = 145 ± 8 mV/s. 

It can be noticed that this law in fact reproduces to a good degree of approximation (error 
smaller than 10%) the two straight lines in the given intervals described, and interpreted, 
above. We can also see that it is verified, as it should, a saturating behavior, namely, 
the infoentropic index a tends to the value 1 for large v, the minimum value it can take. 
Hence, according to Eq.Q it follow that dj tends to 2, what can be interpreted, as already 
noticed, that as w increases the charges keep "feeling" less and less the microroughness at the 
interface. Hence, for v very large there follows that ^ 2 (an averaged smooth surface), 
while the limit of v going to zero would evidence the more detailed effect of the fractal 
topography, given in this case the values df ~ 2.8 (then near 3 and then indicating a strong 
roughness) and a ^ 1.62 ± 0.08, in accordance with the fact that a must remain in the 
interval as given by Eq. (j25|) . and of course 2 < df < 3. 
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IV. CONCLUDING REMARKS 



We have considered in detail a quite interesting case consisting of a situation present in 
certain classes of systems where it becomes difficult for the researcher to apply the standard 
formalism in the conventional, well established, logically and physically sound Boltzmann- 
Gibbs statistics, because of the presence of inhomogeneous fluctuations which are "anoma- 
lous" in some sense and eventually, as in the case here, involving one type or other of 
fractality. In that cases - as already noticed - it may be convenient to depart from con- 
ventional statistics, introducing alternative ones. In other words, the appropriate use of the 
Boltzmann-Gibbs formalism may in some problems be quite complicated to apply or, and 
more fundamental, we may be faced with an inaccessible characterization constraints (infor- 
mation): thus, simplicity or the force of circumstances, respectively, would require setting 
on the use of other measures of informational entropies leading to heterotypical probability 
distributions 2^ accompanied, as noticed, by the use of their associated escort probability 
4l|, and here we found convenient to use a physical statistical mechanics based 



on Renyi s statistics 
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Ml 



The matter has been illustrated in this paper by considering its application in a theo- 
retical study of transport in electrodes in microbatteries, particularly cyclic voltammetry, 
which provides information on the fractal-like microroughened texture of them, providing an 
excellent illustration for the use and validation of such statistical theory. In ideal conditions 
(quite smooth surface) the peak current in the voltammetry experiment is proportional to 
the square root of the sweep-rate velocity, with which the applied fleld is changing. The 
standard theory relates this result to the diffusion of the charges that are going to reach 
the thin fllm electrode. But experimental data show that it is verifled a power law, with a 
power different than 0.5, and this is referred to as resulting from a process of "anomalous" 
diffusion, as noted in the previous sections. As we have noticed, a description in terms of 
a process of diffusion does not apply, once the hydrodynamic motion that proceeds at the 
interphase electrolyte-electrode where is present a morphology with variations (asperity) 
ranging in the nano- to sub-nanometric scale. Hence, a proper description of the hydrody- 
namic motion requires to go well beyond the long-wavelength limit of classical (Onsagerian) 
thermo-hydrodynamics, being required to introduce an extended higher-order (sometimes 
referred-to as mesoscopic) thermo-hydrodynamics, together with the introduction of com- 
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plicated boundary conditions. 

As shown in the Appendix, using Renyi's statistics but ignoring the nano- to sub- 
nanometric-sized roughness of the thin film electrode, it is derived an "anomalous" (non- 
Fickian) diffusion equation, which allows to describe the experimental results in cyclic 
voltammetry. The relevant point is that the index a, that the Renyi statistics contains, 
determines the power law in the voltammetry experiment which allows to estimate the 
roughness of the electrode surface and, then, its quality and influence on the functioning of 
the device. 

In section III has been presented a theoretical analysis of results in a given experiment 
(cyclic voltammetry and atomic force microscopy): in tables and figures are shown the results 
of the study of the index a of Renyi statistics, together with other power indexes (determined 
by the former) and the statistical fractal dimension, all depending on the experimental 
protocol. 

Closing this section we notice that the use of an alternative statistics even though does 
not provide a complete physical picture of the problem in hands - once we do have an 
utter difficulty to handle the appropriate extended non-classical thermohydrodynamics to 
be applied - the sophistication of the formalism allow us to obtain a good insight into the 
physical aspects of the situation. On the one hand, as noticed, the use of the escort prob- 
ability takes care of introducing correlations (fluctuations and higher-order variances), and 
on the other hand the heterotypical probability distribution (the one of Renyi here) modifies 
the weight of the Fourier amplitudes ln(Q)| of the density in relation to their classical val- 
ues in the linear Onsagerian regime (4 71 . l48| . As previously noticed and discussed, classical 
thermo-hydrodynamics together with Fick's diffusion equation is a satisfactory approach 
while |^^(Q)| has leading contributions for small Q (large wavelengths) and negligible for in- 
termediate to large wavenumbers (intermediate to short wavelengths). When not only small 
wavenumbers but ever increasing ones become relevant for the description of the motion, 
i.e. associated to non-negligible amplitudes |^^(Q)| contributing in the Fourier analysis of 
the density, we need to introduce a higher-order thermo-hydrodynamics ( with now equa- 
tions of evolution of the Maxwell- Cattano-type, Burnett and super-Burnett-type, and so 



on 



3, 



4j|. In the case we have analyzed, Renyi heterotypical distribution, which leads to 



obtain the "anomalous" diffusion equation of Eq.(24), and where a > 1, in the restricted 
hydrodynamic description used, decreases what would be the values of |?t-(Q)| in the stan- 
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dard treatment for short wavenumbers, that is, in the domain of classical hydrodynamics, 
while increases those amplitudes outside that region, i.e. in the regime of the higher order 
thermo- hydrodynamics. 
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APPENDIX A: "ANOMALOUS" (NON-FICKIAN) DIFFUSION 

Let us consider first the conventional case of diffusion, when the principle of sufficiency is 
satisfied, meaning that the several stringent restrictions its validity requires are met, namely, 
local equilibrium, linear Onsager relations and symmetry laws, condition of movement with 
long wavelengths and very low frequencies, and weak fluctuations, are verifled. A specific 
criterion of validity is given in Ref. [46], where it is considered a system composed of two 
ideal fluids in interaction between them. The continuity equation for the flux In{r,t) [cf. 
Eq.(19)], after transforming Fourier in time, takes the form: 

{1 + iu Tnl)In{r,t) + TnlV ■ li^\r,t) = , (Al) 

where x„i is the momentum relaxation time. But, at very low frequencies ujTni ^ 1, and 
then 

I„(r,t)^-7:„iV-/(')(r,t), (A2) 

and a direct calculation tell us that 

V./[2l(r,t) = My^(r,t), (A3) 
m 

Replacing Eq. ()A3|) in Eq. ()A2|) . and the latter in the conservation equation for the density, 
(Eq.(14)) we obtain the usual Pick's diffusion equation 

|-n(r,t) + DV'n(r,t) = 0, (A4) 
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where D = ^v^f^ r„i, with \mv'^j^ = ^kBT; Vth is the thermal velocity and D is the diffusion 
coefficient, with dimensions cm^/sec, and m is the mass of each particle. 

Let us now go over the unconventional treatment, which is required once one is looking 
forward for an analysis of data on the basis of a description in terms of a diffusive movement, 
when this is not possible, as a consequence that diffusion in the microroughened region 
is governed by not too long wavelengths (up to nanometric ones, i.e. 10"'' cm, while the 
limitation of the diffusive region is of the order of D/vth, say, typically 10 ^ to 10 ^ 



cm). A higher-order thermo-hydrodynamics 



43 
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needs be introduced, but if the 



lower order description including only the density and its flux is kept, then sufficiency is not 
satisfied and we need to introduce auxiliary Statistical Mechanics [2^ . Let us consider 
the auxiliary statistical operator which for this system of identical free particles is in Renyi 
statistics given by 



where 



(A5) 



Ani(r, p; t) = ni(r, p; |r) - (ni(r, p; t)) 



a ' 



(A6) 



with 



(fii(r,p; t))^ 



rfrfii(r,p;|r) V^{T\t) 



(AT) 



and where (/9a(r, p,t) is the associated Lagrange multiplier, rj ensures its normalization, fii 
is given in Eq.(9), and 



dV {q^)' 



(A8) 



is the accompanying escort probability (cf. Eq.(l)). 
Introducing the modified Lagrange multiplier 



V?a(r,p,t) = v?a(r,p,t) 



l-{a-l) d-'r d'p ifair, p, t) (ni(r, p; t))^ 



-1 



(A9) 



we find that 



l + (a-l) / d^r / rf>(^„(r,p,t)ni(r,p|r) 



(AlO) 
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with z{t) ensuring the normahzation condition. At this point we introduce for (p the form 

2 

(^„(r, p, t) = Fh{r, t)^ + F„(r, t) + F„(r, t) ■ ^ , (All) 
Im m 

for, in a such a way, keeping as basic variables the three of Eqs.(lO) to (12), and thus arriving 
to the statistical operator of Eq.(13). 

Using Eqs. (A. 5) to (A. 14), after some lengthy but straightforward calculations, it follows 
for the energy density that 

h{v,t) = Jd^P^J rfrni(r,p|r)P,(r,t) = M(r,t) + n(r,t)^mt;2(r,t), (A12) 

i.e., composed of the energy associated to the drift movement (the last term) and the internal 
energy density 

u{r,t) = -^^^n^-{r,t), (A13) 
5 -3a p^{r,t) 



where 

2 Q-l 

r 9-Tr N r 9 TTJ 1 9. n, 9, ^ 

Ca{r,t) 



(a- 1)3/2^1 (t) [4(r,t) 



• (A14) 



i3(z/, x) is the Beta function, we have written F„Q,(r, t) = /?Q,(r, t); F„Q,(r, t) = m /^^(r, t) v(r, t) 
(introducing a "drift velocity" field v(r, t)); moreover 

5 — 3 a / A . 

7a = 7^ F ' A15 

3 a — 5 

and the values of a are restricted to the interval 

1 < « < ^ • (A16) 



Finally, the second order flux is given by 

/[fl(r,t) = jd'P^^I dTMr,p\T)V^{T,t) 



= n(r, t), [v,(r, t) v,(r, t)] + ^(r, t) l^l , (A17) 

3m 

where l'^' is the unit second- rank tensor, [...] stands for the tensorial product of vectors, 
and it can be noticed that 

pPl(r,t) = m/P](r,t)-mn(r,t) [v,(r, t) v„(r, t)] , (A18) 
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is the pressure tensor. Neglecting the terms quadratic in the drift velocity, combining the 
above equations we obtain that 

/Pl(r,t)=U(r,t)n^'^(r,t)lP], (A19) 

where 

is the quantity present in Eq.(22). 



2 5 Ca{r,t) 
3m 5 - 3a ^^^(r,^) 



(A20) 



APPENDIX B: FLUX RELAXATION 



The scattering integral of Eq.(16) depends, through the statistical operator V^, on the 
non-equilibrium thermodynamic quantities Fh, Fn, and F„ (related to the Lagrange multi- 
pliers that the variational method introduces; cf. Appendix A). But the latter are related 
to the basic variables, h, n and I„, by means of Eqs.(lO) to (12), which are considered 
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as being equations of state in the associated nonequilibrium thermodynamics 
Hence J„Q,(r, t) is a functional of these basic variables. But we do have, on the one hand, 
that n(r, t) = no + A?2(r, t) and /i(r, t) = /iq + A/i(r, t), namely their constant global values 
modified by a small variation, that is An ^ uq and, on the other hand, I being also small 
we can write a series expansion in it but retaining only the first (linear) contribution 

J„,(r,t)~ [eL2]]-'®I„(r,t), (Bl) 

where [0a ^] ^ is the second-rank tensor (playing the role of the inverse of a tensorial 
Maxwell-characteristic time j^) of components 5Ji/5Ij, dependent only on riQ and h^. 
Writing 

'la 

where 

1 1 



and the last term on the right of Eq.(B2) being then the traceless part of Maxwell charac- 
teristic time tensor, neglecting it Eq.(Bl) becomes Eq.(18). 
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Cyclic Voltammetry 





e 


a 




slow scan rates 


0.85 


1.58 


2.70 


fast scan rates 


0.63 


1.34 


2.26 


Atomic Force Microscopy 






r 


df 


short L ranges 




0.76 


2.24 


medium L ranges 




0.27 


2.73 



TABLE I: Exponent e, Infoentropic index a and fractal dimension df obtained from cyclic voltam- 
metry experiments performed at slow and fast scan rates. The values of the roughness exponent r 
and fractal dimension df obtained from Atomic Force Microscopy images at short and medium L 
ranges are also presented. 
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Figure Captions 

Figure 1- Log-log plot of peak current ip versus sweep rate v. 

Figure 2- AFM micrograph, 1000 nm x 1000 nm x 23 nm. 

Figure 3 - Log-log plot of interface width W{L) versus partitions L. 
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FIG. 1: Log-log plot of peak current ip versus sweep rate v. 
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FIG. 3: Log-log plot of interface width W(L) versus partitions L. 
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